Methods to print N-tone images with multi-leveling techniques

ABSTRACT

Generating a N-tone image from a grey scale image with a number of multi-level dither matrices. Both images have many pixels, with each pixel in the grey scale image having a grey level that is within a maximum and a minimum levels, and with each pixel in the N-tone image having a level that is one of N levels with N bigger than two. The method includes the step of determining the level of each pixel in the N-tone image based on the levels of its corresponding pixel in the grey scale image, and an original dither matrix, the separations between adjacent levels in the N levels are not uniform. After the determination process, every pixel in the N-tone image is modelled to represent its level. Based on the modelling, the N-tone image is printed.

CROSS REFERENCE TO RELATED APPLICATION

This is a continuation of application Ser. No. 08/205,672 filed on Mar. 2, 1994, now abandoned.

CROSS REFERENCE TO RELATED APPLICATION

The present invention is related to commonly assigned and co-pending U.S. application entitled, "Print N-Tone Images With Imperceptible Dots", filed Mar. 2, 1994, and invented by Qian Lin.

BACKGROUND OF THE INVENTION

The present invention relates generally to printing an image and more particularly to printing an image with N-tones.

Numerous methods are available to print a halftone image from a grey scale image. The methods usually involve establishing the approximate grey level of each pixel of the grey scale image, and then, based on some representation schemes, printing dots to represent the grey scale image.

One form of representation scheme depends on a dither matrix, which has the same number of pixels as the grey scale image. Each pixel in the matrix has a level, which is compared to the level of its corresponding pixel in the grey scale image to produce the level of a pixel in the halftone image. A general discussion of a dither matrix to render an image can be found in "Digital Halftoning," by R. Ulichney (1987). Another form of representation scheme depends on the error diffusion technique, with a general discussion found in "An Adaptive Algorithm for Spatial Greyscale," written by Floyd and Steinberg, and published in the Proc. SID, Volume 17, pages 75-77, 1976.

One type of dither matrix is known as a Bayer matrix. FIG. 1 shows a prior art grey ramp at 300 dots per inch with 256 levels, printed using a Bayer matrix with 8 by 8 pixels. More than 50% of the grey levels in the ramp are indistinguishably black. This reduces the number of grey levels to print an image. Also, many transitions between distinguishable grey levels are not smooth, causing conspicuous discontinuities in a figure generated by the Bayer matrix.

It should be apparent from the foregoing that there is still a need for a method to print, relatively inexpensively, an image with more distinguishable grey levels and smoother transitions between adjacent levels.

SUMMARY OF THE INVENTION

The present invention is on methods of generating, from a grey scale image, a N-tone image with more distinguishable grey levels and better transitions between adjacent levels, in a relatively inexpensive manner. The method is based on generating a N-tone image, instead of a halftone image.

The grey scale image and the N-tone image have many pixels. Each pixel in the grey scale image has a level that is within a maximum and a minimum level. Each pixel in the N-tone image has a level that is one of N levels, with N bigger than two.

In one preferred embodiment, the level of each pixel in the N-tone image is determined based on the levels of its corresponding pixel in the grey scale image and based on an original dither matrix; the separations between adjacent levels in the N levels are non-uniform. After the determination process, every pixel in the N-tone image is modelled to represent its level. Based on the modelling, the N-tone image is printed.

One preferred method to model each pixel is to replace every pixel in the N-tone image by many contiguous sub-pixels. The level of each pixel is represented by the intensity for each of its corresponding sub-pixels so that different combination of intensities of the sub-pixels produce different levels for that pixel. Based on the above method, the N-tone image printed has smoother transitions between adjacent grey levels and more distinguishable grey levels in its grey ramp.

Other aspects and advantages of the present invention will become apparent from the following detailed description, which, when taken in conjunction with the accompanying drawings, illustrates by way of example the principles of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a prior art grey ramp printed using a Bayer matrix.

FIG. 2 shows a representation of the first preferred embodiment of the invention.

FIG. 3 shows a preferred table of groups of grey levels to generate a N-tone image for the first preferred embodiment.

FIG. 4 shows one preferred way to form a multi-level dither matrix in the present invention.

FIG. 5 shows one multi-level dither matrix formed by the method described in FIG. 4.

FIG. 6 shows one preferred way to form the N-tone image of the first preferred embodiment.

FIG. 7 shows one preferred way to generate subpixels and segments in the first preferred embodiment.

FIG. 8 shows the formation of sub-pixels in the first preferred embodiment.

FIG. 9 shows one preferred sets of segments for four sub-pixels.

FIG. 10 shows the 16 levels printed using the different segments in FIG. 9.

FIG. 11 shows the relationship between the grey levels and the density of each level in the N-tone image of the first preferred embodiment.

FIG. 12 shows a picture printed by the first preferred embodiment.

FIG. 13 shows an example for the second preferred embodiment, with 16 contiguous sub-pixels to form one pixel.

FIG. 14 shows the turn-on sequence of the subpixels to generate 13 different levels for the second preferred embodiment.

FIG. 15 shows one table of groups of grey levels for the second preferred embodiment.

FIG. 16 shows the grey ramp generated based on the second preferred embodiment of the present invention.

FIG. 17 shows a part of a fourth preferred embodiment of the present invention.

FIG. 18 shows a preferred way for the fourth preferred embodiment.

FIG. 19 shows graphically one preferred way to form the special dither matrix in the fourth preferred embodiment.

FIG. 20 shows a set of preferred steps to form the special matrix in FIG. 19.

Same numerals in FIGS. 1-20 are assigned to similar elements in all the figures. Embodiments of the invention are discussed below with reference to FIGS. 1-20. However, those skilled in the art will readily appreciate that the detailed description given herein with respect to these figures is for explanatory purposes as the invention extends beyond these limited embodiments.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 2 shows a representation of the first preferred embodiment 100 of the invention, which prints a N-tone image 102 by a printer, a plotter or other imprinting device 104 from a grey scale image 106. The N-tone image 102 has two intermediate forms, the first intermediate N-tone image 101 and the second intermediate N-tone image 103. The images have many pixels; for example, the grey scale image 106 has pixels 108 and 112; and the first intermediate N-tone image 101 has pixels 114 and 116.

Each pixel in the grey scale image has a level that is within a maximum and a minimum level. For example, the maximum level is 255 and the minimum level is 0; each level represents one grey level of the grey scale image 106.

The first intermediate N-tone image 101 has N levels, with N preferably greater than 2. The level of each pixel in the first intermediate N-tone image 101 is from one of the N levels. The N levels are selected from the grey levels in the grey scale image 102 by dividing the grey levels into N groups. Each group of grey levels has a group maximum level, which is the maximum level in the group and is one of the N levels. The separations between adjacent levels in the N levels are non-uniform. The method to set the non-uniformity will be described later in the specification.

FIG. 3 shows an example of a preferred table of groups of grey levels to generate a N-tone image with 16 levels. The sixteenth group 202 covers the range of grey levels between 255 to 253 with its group maximum level being 255. Every group is represented by its group maximum level. The entire 16 groups give rise to 16 group maximum levels, which represent the 16 levels of the N-tone image.

The first preferred embodiment includes a number of multi-level dither matrices 118, such as the multi-level dither matrix 122. Each multi-level dither matrix has many pixels, such as the matrix 122 has pixels 124 and 126.

The multi-level dither matrices are formed through the groups of levels shown in FIG. 3 and an original dither matrix. FIG. 4 shows one preferred way 260 to form a multi-level dither matrix 122 as the one shown in FIG. 5. First, N groups of levels are formed, 262, as in FIG. 3. Then the matrices are formed through scaling 264 the original dither matrix according to the groups of levels. The original dither matrix 250 may be a Bayer matrix or other types of dither matrix with levels ranging such as from 1 to 255; the matrix is the type that has been used extensively to render an image and will not be further described in this application. The multi-level dither matrix 122 shown in FIG. 5 corresponds to the fourth group 206 of grey levels shown in FIG. 3, which covers the grey levels from 56 to 31. To form the multi-level dither matrix 122, the levels in the original dither matrix (1 to 255) are scaled according to the levels in the group 206, with grey levels ranging from 56 to 31; for example, the level 167 in the second pixel of the original dither matrix becomes the level 47 in the second pixel of the multi-level dither matrix using the following direct scaling calculation, with fractions rounded off:

(((56-31)/(255-1))*167 +31).

Similarly, the level 40 in the pixel 257 of the original dither matrix becomes 35 in the pixel 259 of the multi-level dither matrix. Based on the above method, the multi-level dither matrix 122 is formed. Using the original dither matrix and the 16 groups of grey levels shown in FIG. 3, sixteen multi-level dither matrices are generated by the direct scaling calculations. FIG. 6 shows a set of preferred steps to form the N-tone image, based on the multi-level dither matrices.

To generate the first intermediate N-tone image 101, every pixel in the grey scale image 106 is compared to and thresholded by a pixel in one of the multi-level dither matrices. For example, the first pixel 108 in the grey scale image 106 has the grey level 51. This grey level falls within the fourth group of grey levels shown in FIG. 3. The multi-level dither matrix 122, generated by the fourth group of grey level, is selected, 275. The first pixel of the grey scale image is then compared, 277, to the first pixel of the multi-level dither matrix 122. In the comparison process, first, one decides if the level of the pixel in the grey scale image is larger than or equal to the level in the corresponding pixel of the multi-level dither matrix; if it is, the group maximum level of the corresponding group generating the multi-level dither matrix is chosen for the first intermediate N-tone image; however, if the level of the pixel in the grey scale image is smaller, the group maximum level of the subsequent group is chosen for the first intermediate N-tone image. In the present example, with 51 bigger than 50, the group maximum level, 56, is chosen to be the level of the first pixel in the first intermediate N-tone image. For the second pixel, with 46 less than 47, the group maximum level of the subsequent group, 31, is chosen to be the level in the N-tone image. From the levels in the grey scale image and the multi-level dither matrices, the pixel-to-pixel comparison method, as described above, generates the first intermediate N-tone image 101.

The second intermediate N-tone image 103 is generated by two steps. First, each pixel in the first intermediate N-tone image 101 is represented or modelled, 279, by a number of contiguous sub-pixels, and then each sub-pixel is represented or modelled by a number of segments.

While not wishing to be bounded by theory, it is believed that some of the advantageous results of the invention are obtained through understanding the functionality of relatively low-cost printers. To generate an image with substantially imperceptible dots, every pixel is printed by, for example, a laser printer with a resolution of at least 600 pixels or dots-per-inch. The characteristics of each pixel is preferably controlled through a pulse with a pulse width. The pulse width can be sub-divided into a number of segments. Each segment can be at a high or a low intensity; at a high intensity means that the segment is dark or that the segment is turned on with the pulse width extending into the segment; and at a low intensity means the segment is light or the segment is turned off with the pulse width not extending into the segment. For a typical present-day 600 dots-per-inch laser printer, if the pulse is divided into eight segments, it is preferable to have more than one segment at a high intensity to generate a reproducible output. In other words, every segment at a high intensity preferably should have a neighboring segment also at a high intensity. Based on this representation, different levels of the image require different pulse widths or require selecting different number of segments. In order to generate more levels for the output image, every pixel in the N-tone image is represented or modelled by a number of pixels; the number of pixels is designated as sub-pixels. In one preferred embodiment, by combining four sub-pixels to generate different levels for each pixel, one gets N-tone images with a 300 dots-per-inch resolution; the N-tone images generated have substantially imperceptible dots for the unaided eyes of an ordinary person, with 20/20 vision, under normal condition. The above theory also applies to other relatively low-cost printers, such as ink-jet printers, by varying the intensity of each printed dot in a way similar to varying the pulse width in each pixel of a laser printer. Then, using the invented methods in the present invention, one may be able to generate images with substantially imperceptible dots.

One example of using sub-pixels and segments is illustrated in FIGS. 8 and 9, with a set of preferred steps shown in FIG. 7. The example is based on a 600 dots-per-inch laser printer. FIG. 8 shows each pixel in the first intermediate N-tone image 101 being replaced, 283, by four contiguous sub-pixels; for example, the pixel 114 is replaced by the sub-pixels 302, 304, 306 and 308. FIG. 9 shows the four sub-pixels, with each sub-pixel replaced, 285, by eight segments. There is a number within each segment; for example, the segment 310 has a number of 15, and the segment 312 has a number of 1. The number inside the segment represents the level of the pixel when the segment is turned on or is at the high intensity. For example, for the level 15, the pulse width extends across three segments, the segment 310, 314 and 316; for the level 1, the darkest level, the pulse width extends across all the segments; and for the level 16, no segment is turned on or all segments are at the low intensity. Thus, the level of each pixel is represented, 287, by the segments and the sub-pixels.

The sixteen levels have perceptibly different levels of grey. One preferred way to decide which segment to be turned on for which level is based on one's visual perception. In the present embodiment shown in FIG. 9 with 32 segments, the segments are turned on one after the other, with the level created measured by a densitometer, which measures the density of a level. Theoretically, there could be 33 levels from the 32 segments, but many levels are visually indistinguishable. Those indistinguishable levels are grouped as one level, for example the numerous segments grouped together for level 1 and level 2. In this embodiment, sixteen distinguishable levels are selected. Another observation in the preferred embodiment, using a 600 sub-pixels-per-inch laser printer, is that the lightest level 15 preferably needs more than 1 segments. The effect of a segment turned on by itself without any neighboring segments being turned on may not be very reproducible. Thus, the lightest level preferably needs more than one segment to be turned on. However, once there are segments turned on, subsequent levels may only need one additional segment to be perceptibly different. For example, the level 14 only has one more segment turned on, as compared to the level 15. In this embodiment, the number of segments between levels are not uniform; for example, the difference between levels 16 and 15 is 3 segments, and the difference between levels 15 and 14 is 1 segment only.

Another preferred way to decide which segment to be turned on for which level is to retain the 33 levels for the 32 segments. A number of levels may not be visually distinguishable, and one needs 6 bits to designate all the levels. If the number of bits is of concern, one may remove one level, and then only 5 bits would be needed to designate the 32 levels.

The above example is based on a 600 sub-pixels or dots-per-inch laser printer. Laser printers with higher resolution can be used. The number of sub-pixels for each pixel does not have to be 4. The number of segments for each sub-pixel does not have to be 8. By experimenting with the number of segments, the number of sub-pixels and the dot size, one, based on the teachings in the present disclosure, can generate images with substantially imperceptible dots, with a different resolution, a different number of sub-pixels and a different number of segments.

FIG. 10 shows the 16 levels printed out using the different pulse widths or segments shown in FIG. 9. For example, the fourth square block, 325, is the fourth level. As shown in FIG. 9, all the segments having number 4 or higher in it will be turned on, which means that pulses are turned on for three sub-pixels, 304, 304 and 306; in other words, 24 of the 32 segments are turned on.

The 16 levels shown in FIG. 10 are measured by the densitometer. The density of each level of the N-tone image is related back to the grey levels of the grey scale image. In one preferred embodiment, the relationship is assumed to be linear. FIG. 11 shows the relationship between the grey levels and the density of each level of the N-tone image. Each level in FIG. 10 is mapped to a group or a range of grey levels. For example, the fourth level approximately has a grey level of 56, and the third level has 31; so the fourth group of levels covers the grey levels from 56 to 31. This is one preferred method to determine all the groups of grey levels in FIG. 3.

From the second intermediate N-tone image 103, the N-tone image 102 is printed, 281, by the printer 104. In the present embodiment, one way to represent a segment being turned on is to use a value of 1, and a segment being turned off by a value of 0. The segments with 1 in them will be printed. FIG. 12 shows a picture printed out by the first preferred embodiment. The image has 300 dots or pixels per inch. Each pixel is represented by 4 sub-pixels, and each sub-pixel by 8 segments. There are altogether 32 segments for each pixel. In this embodiment, each pixel only has 16 levels, which can be represented by 4 bits. Thus, with an increase of 4 bits for every pixel, for an ordinary person, the figure generated has substantially imperceptible dots. Any dots in the figure, such as those on the apple, are in the original grey scale image.

A second preferred embodiment of the present invention is similar to the first except one does not have to vary the pulse width, or use different segments for each sub-pixel in the second intermediate dither matrix 103. The N-tone image is generated, for example, through an original Bayer matrix with 8 by 8 pixels.

In one example of the second embodiment, every pixel is represented or modelled by 16 contiguous sub-pixels, as shown in FIG. 13. FIG. 14 shows the turn-on sequence of the sub-pixels to generate 13 different levels. For example, if one prefers the fourth level, every sub-pixel labelled 4 or higher will be turned on, or will have a dot printed in it. This turn-on sequence is determined based on a classical screen or cluster dot at 45° model; such model should be obvious to those skilled in the art and will not be further described here. From the 13 levels, a densitometer generates 13 groups of levels according to a graph similar to FIG. 11. FIG. 15 shows the groups of levels found. Again, the separations between the groups of levels are not uniform.

From the 13 groups of levels in FIG. 15, the original Bayer matrix forms 13 multi-level dither matrices. From the multi-level dither matrices, the first intermediate dither matrix 101 is generated. Based on the 13 levels as shown in FIG. 14, the second intermediate N-tone image and the N-tone image are produced as in the first preferred embodiment.

FIG. 16 shows the grey ramp with 256 grey levels generated based on second preferred embodiment with the multi-level Bayer matrices and the 13 levels. The pixels have a resolution of 75 dots-per-inch, with the subpixels at 300 dots-per-inch. As compared to the prior art grey ramp in FIG. 1, there are significantly more distinguishable grey levels. Moreover, the transition between grey levels in FIG. 16 are more smooth than those in FIG. 1. These effects are more obvious to a viewer if the two figures are viewed side-by-side at a distance, such as 2 feet away from the viewer.

For a third preferred embodiment, its main difference from the first preferred embodiment is that the multi-level dither matrices are replaced by a multi-levelling error diffusion technique to generate the first intermediate N-tone image 101. In the normal error diffusion technique, the error from each pixel is compared to the middle grey scale level of the grey scale image, with the errors diffused to its surrounding pixels. Error diffusion techniques are well-known to those with ordinary skill in the art and will not be further described in the specification. For the preferred multi-level error diffusion technique in the present invention, the level in each grey scale image pixel again maps to a corresponding group of grey levels, such as the one shown in FIG. 3. The level in each pixel of the grey scale image is compared to the middle level of its corresponding group of levels, and the error is again diffused to its corresponding pixels to generate the first intermediate N-tone image 101. For example, the first pixel 108 with a level of 51 is compared to the level 43 (the average of 56 and 31), and the error is diffused to its neighbors.

FIG. 17 shows a portion of a fourth preferred embodiment 401 of the present invention, with a set of preferred steps shown in FIG. 18. The difference between the first and the fourth embodiment is that the fourth embodiment reduces the number of steps in the first embodiment by collapsing many steps into one special dither matrix. In this embodiment, the level of each segment in a second intermediate N-tone image 404 is determined by a pixel-to-pixel comparison between the grey scale image 400 and the special dither matrix 402. Based on the comparisons, the level in each pixel of the N-tone image is determined, 375. Then, from the second intermediate N-tone image 404, the N-tone image is printed, 377, by the printer. In another embodiment, one preferably does not need to determine, 375, the level in each pixel. All the segments in the second intermediate N-tone image 404 are sent as a bit map to the printer, and the bit map is printed directly.

In one preferred embodiment, as shown in FIG. 17, every pixel in the grey scale image 400 is subdivided into a group of contiguous segments to form a finer-resolution grey scale image 406. All the segments for each pixel may have the same grey level as its corresponding pixel in the grey scale image 400; for example, the first pixel 408 with a level of 206 is subdivided into thirty-two contiguous segments, such as 410, 412, 414 and 416; all the segments have the level of 206. For another embodiment, the finer-resolution grey scale image 406 actually has a finer resolution than the grey scale image 400; for example, thirty-two times higher in resolution. After the sub-division, each segment in the finer-resolution grey scale image is compared to its corresponding segment in the special dither matrix 402.

FIG. 19 graphically shows one preferred way to form the special dither matrix 402, with a set of preferred steps shown in FIG. 20. As an overview, many standard grey scale images, such as 460, 462 and 464, are first compared, 500, to the multi-level dither matrices, such as 118 and 122, to produce many intermediate output matrices, such as 450, 452 and 454. Then, every pixel in the intermediate output matrices is replaced by a group of contiguous sub-pixels as in the first preferred embodiment to form, 502, many output matrices, such as 470, 472 and 474. All the output matrices are added, 504, together to form the special dither matrix 402.

In more detail, every pixel in each standard grey scale image has the same grey level. For example, all the pixels in the standard grey scale image 460 has the grey level 0. The grey scale image may have 256 levels. In order to cover all the levels, there are preferably 256 standard grey scale images.

The multi-level dither matrices are generated as in the first preferred embodiment. For each standard grey scale image, its pixels are compared to one corresponding multi-level dither matrix to generate its intermediate output matrix. The comparison is similar to the comparison in the first preferred embodiment. For example, the standard grey scale image 462 has the level of 45 for all its pixels. Its corresponding multi-level dither matrix is 122, which covers the levels from 56 to 31. Thus, every pixel in the standard grey scale image 462 is compared to its corresponding pixel in the multi-level dither matrix 122 to generate the intermediate output matrix 452.

All the intermediate output matrices are transformed to their corresponding output matrices by replacing every pixel with a group of contiguous segments. This replacement process is similar to the replacement process in generating the second intermediate N-tone image from the first intermediate N-tone image in the first embodiment. For example, each pixel in the intermediate output matrices is replaced by 32 segments with their corresponding intensity to represent the level of each pixel; the segments that should be turned on have one in them, and those that should be turned off have zero. After the replacement, every intermediate output matrix becomes its corresponding output matrix. For example, the intermediate output matrix 452 becomes the output matrix 472 after the replacement process. All the output matrices are then added together by matrix addition to produce the special dither matrix 402.

A Bayer matrix may be used for the fourth preferred embodiment. The difference is that every pixel in the intermediate output matrix is not replaced by segments, but by sub-pixels only, as in the second preferred embodiment. A sub-pixel that should be turned on has one in it and the one that should be off has a zero in it.

From the foregoing it should be appreciated that the present invention generates visually pleasing N-tone images in a relatively inexpensive manner. The above description can be extended to color images. For a color N-tone image, the invented method is repeated at least two more times, each time for a different color to generate the color N-tone image. The present invention can also be extended to a display with the N-tone image shown on the display, instead of printed on a printer. In fact, the printing step in the present invention includes displaying, with the printer in the present invention being a display.

Other embodiments of the invention will be apparent to the skilled in the art from a consideration of this specification or practice of the invention disclosed herein. It is intended that the specification and examples be considered as exemplary only, with the true scope and spirit of the invention being indicated by the following claims. 

I claim:
 1. A method of generating N-intensity-tone image from a grey scale image, both images having a plurality of pixels, each pixel in the grey scale image having a grey level that is within a maximum and a minimum levels, and each pixel in the N-intensity-tone image having a grey level that is one of N intensity levels with N bigger than two, the method comprising the steps of:determining the intensity level of each pixel in the N-intensity-tone image based on the level of its corresponding pixel in the grey scale image and based on an original dither matrix, the number of grey levels between adjacent levels in the N intensity levels being not constant; modeling the intensity level of every pixel in the N-intensity-tone image; and printing the N-intensity-tone image with a printer by printing the intensity level of each pixel based on the modeling; such that: the printer prints the intensity level of each pixel by controlling a pulse with a pulse width; and the pulse width of one intensity level is different from the pulse width of at least one other intensity level.
 2. A method as recited in claim 1 wherein the step of determining further comprises:selecting for every pixel of the grey scale image a multi-level dither matrix from a plurality of multi-level dither matrices, which are generated based on the N-intensity levels and the original dither matrix; and comparing the level of each pixel in the grey scale image and the level of the corresponding pixel in the selected multi-level dither matrix to determine the level of the corresponding pixel in the N-intensity-tone image.
 3. A method as recited in claim 2 wherein there are N multi-level dither matrices, which are generated by the steps of:dividing the levels in the grey scale image into N groups based on the N intensity levels of the N-intensity-tone image; and scaling the original dither matrix by the levels in each group to generate the N multi-level dither matrices.
 4. A method as recited in claim 1, wherein there are visually perceptible change in intensity between adjacent levels in the N intensity levels of the N-intensity-tone image.
 5. A method as recited in claim 1, wherein the N-intensity-tone image is in black-and-white.
 6. A method as recited in claim 1, wherein:the N-intensity-tone image is in color; and the method further comprises the step of repeating the method more than once, each time for a different color to generate the color N-intensity-tone image.
 7. A method as recited in claim 3, wherein:the N-intensity-tone image is in color; and the method further comprises the step of repeating the method more than once, each time for a different color to generate the color N-intensity-tone image.
 8. A method as recited in claim 1 wherein the step of modelling includes:replacing every pixel in the N-intensity-tone image by a plurality of sub-pixels, each sub-pixel having an intensity; and representing the level of each pixel by the intensity of each of its corresponding sub-pixels.
 9. A method as recited in claim 1 wherein the step of printing includes using a black-and-white printer.
 10. A method as recited in claim 1 wherein the step of printing includes using a color printer.
 11. A method of generating N-intensity-tone image from a grey scale image based on an original dither matrix and a plurality of standard matrices, the matrices and the images having a plurality of pixels, each pixel in the grey scale image having a grey Ievel that is within a maximum and a minimum levels, and each pixel in the N-intensity-tone image having one of N levels with N bigger than two, the number of grey levels between adjacent levels in the N intensity levels being not constant, the method comprising the steps of:determining the intensity level of each pixel in the N-intensity-tone image by comparing the level of its corresponding pixel in the grey scale image to the corresponding pixel in a special dither matrix, the special dither matrix being dependent on (a) the N intensity levels, (b) the original dither matrix and (c) modeling with sub-pixels to generate N intensity levels; and printing the N-intensity-tone image with a printer by printing the intensity level of each pixel according to its determined intensity level; such that: the printer prints the intensity level of each pixel by controlling a pulse with a pulse width; and the pulse width of one intensity level is different from the pulse width of at least one other intensity level.
 12. A method as recited in claim 11, wherein the special dither matrix is generated by:operating each standard matrix by one of a plurality of multi-level dither matrices to generate its corresponding intermediate output matrix, the multi-level dither matrices being formed from the N intensity levels and the original dither matrix; producing a plurality of output matrices through modelling every pixel in the intermediate output matrices by the plurality of contiguous sub-pixels; and adding the output matrices.
 13. A method as recited in claim 12, wherein:the pixels in each standard matrix have the same grey level, and the levels from all the standard matrices cover all the levels in the grey scale image; and the multi-level dither matrices are generated bydividing the levels in the grey scale image into N groups based on the N intensity levels of the N-intensity-tone image, and scaling the original dither matrix by the levels in each group to generate the N multi-level dither matrices.
 14. A method as recited in claim 1, wherein:the step of modelling comprises modelling with sub-pixels to generate N intensity levels; and the steps of determining and modelling depend on a special dither matrix and a plurality of standard matrices.
 15. A method as recited in claim 14, wherein the special dither matrix is generated by:operating each standard matrix by one of a plurality of multi-level dither matrices to generate its corresponding intermediate output matrix, the multi-level dither matrices being formed from the N intensity levels and the original dither matrix; producing a plurality of output matrices through modelling every pixel in the intermediate output matrices by the plurality of contiguous sub-pixels; and adding the output matrices.
 16. A method as recited in claim 15, wherein:the pixels in each standard matrix have the same grey level, and the levels from all the standard matrices cover all the levels in the grey scale image; and the multi-level dither matrices are generated bydividing the levels in the grey scale image into N groups based on the N intensity levels of the N-intensity-tone image, and scaling the original dither matrix by the levels in each group to generate the N multi-level dither matrices. 